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Abstract: We show how to generalize the classical duals found by Gabadadze et al to 
a very large class of self- interacting theories. This enables one to adopt a perturbative 
description beyond the scale at which classical perturbation theory breaks down in the 
original theory. This is particularly relevant if we want to test modified gravity scenarios 
that exhibit Vainshtein screening on solar system scales. We recognise the duals as be- 
ing related to the Legendre transform of the original Lagrangian, and present a practical 
method for finding the dual in general; our methods can also be applied to self-interacting 
theories with a hierarchy of strong coupling scales, and with multiple fields. We find the 
classical dual of the full quintic galileon theory as an example. 
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1 Introduction 

For any sensible theory of gravity, the dynamics is described by a complicated coupled 
system of non-linear partial differential equations. This is true of General Relativity (GR), 
and must therefore be true of any theory that hopes to mimic GR at some appropriate scale. 
The complexity of the governing equations renders it difficult to find exact solutions, and 
typically one can only make progress by imposing certain symmetries in order to reduce 
the phase space, or else to consider perturbations about a known background solution. 
Perturbation theory works best far away from the relevant excitation of the background, 
and will inevitably break down as we move closer and closer to the source. In classical GR 
this breakdown occurs at the Schwarzschild radius of the source. This is good enough for 
the most part, since we cannot experimentally probe the dynamics within the Schwarzschild 
radius since it would lie behind an event horizon. 

In modified theories of gravity the situation can be more subtle. Large distance modifi- 
cations of GR are often considered in order to address an outstanding problem in cosmology 
such as the dark matter problem [1], the dark energy problem [2], and/or the cosmological 
constant problem [3] (see [4] for a general review and a more complete list of references). 
This requires an 0(1) deviation from GR on astrophysical or cosmological scales, and for 
the theory to be phenomenologically viable, this deviation must reduce to < 0(1O -5 ) on 
solar system scales. Such dramatic suppression can sometimes be achieved through non- 
linearities, either from non-linear matter couplings, as in the chameleon mechanism [5], or 
via non-linear self-interactions, as in the Vainshtein mechanism [6, 7]; here we are primarily 
interested in the latter. 

Such self couplings have two important effects. At the quantum level, the interactions 
become strongly coupled above some particular scale, A, and one can no longer trust the 
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classical background solution on distances < A -1 . At the classical level, the interactions 
lead to the breakdown of the standard linearized theory around a heavy source. In a 
generic modified gravity theory, with a non-relativistic source of mass M, the linearized 
perturbations break down at the Vainshtein scale, ry, which typically takes the form [8] 



For the Sun, this scale must lie beyond the edges of the solar system, making it much larger 
than the Schwarzschild radius. This is both a blessing and a curse. If the Vainshtein mech- 
anism is effective the non-linearities help to suppress the modifications of GR making the 
theory compatible with observation. However, in the absence of a perturbative description 
below the Vainshtein scale, it is difficult to test any corrections to the leading order effect. 

Recently, Gabadadze et al [9] considered two examples of a classical theory with deriva- 
tive self-interactions, each admitting a standard perturbative description above the relevant 
Vainshtein scale. They then presented a "classical dual" of each theory, describing exactly 
the same physics, but admitting a perturbative description below the Vainshtein scale. This 
is not a duality in the usual sense of strong versus weak coupling. Rather it is a classical 
analogue in the sense that the classical expansion parameter, r/ry, is inverted, with the 
two descriptions giving us the flexibility to do perturbation theory over a much larger range 
of scales, provided the classical effective theory is valid. This is reminiscent of Vainshtein's 
original approach to massive gravity [6] in which he made use of an expansion in ry/r 



graviton mass and r s is the Schwarzschild radius of the source. Indeed, this inversion of 
the expansion parameter was a big factor in motivating Gabadadze et aFs recent work. 

In this paper we give a general prescription for finding the classical duals of a much 
larger class of theories with self-interactions. Our methods work well for those theories 
that remain weakly coupled at low energies, and whose classical high energy dynamics is 
dominated by Appoint interactions with finite A" only. We also require that the interactions 
become subdominant as the fields tend to zero. 

We begin by observing that the classical duals presented in [9] actually make use of 
Legendre transforms of the interaction terms. Running with this idea, we are able to 
generalize their method. Indeed, it is no coincidence that some mathematicians refer to 
the Legendre transform as the Legendre dual[10]. The duality is only useful if it admits a 
new perturbative description, and this is always true for the broad class of theories under 
consideration. 

Whilst the Legendre transform picture is certainly instructive, it is not always straight- 
forward to arrive at a workable dual theory. This is down to technical difficulties in inverting 
the expression for the transformed variables. To alleviate this problem we will also present 
a more practical method for finding the dual, using Lagrange multipliers, which are then 
integrated out. This should really be considered as the working method for finding the 
classical dual of your favourite theory. It also helps us to identify exactly when a useful 
dual can be found. Both methods can be applied to theories with multiple scales, and 
multiple fields, in contrast to the examples given in [9]. 




< a < 1 



i 



above the scale ry ~ 5 > and an expansion in r/ry below the scale ry, where m is the 
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The rest of this paper is organised as follows: We start with a recap of the cubic 
galileon case covered in [9], showing that the dual theory is just the Legendre transform 
of the original, and we then extend this idea in section 3 to the general well as 

proposing a practical implementation in section 3.2. We see how models with multiple 
scales are dealt with in section 3.3 and multiple fields are examined in section 3.4. As 
an example of the method we work out the dual of the full galileon model in section 4, 
allowing for multiple scales, and then finish the paper with a discussion. 



2 The cubic galileon and its classical dual 

Let us begin by reviewing one of the examples considered by Gabadadze et al [9], namely, 
the cubic galileon theory [11], 

C = - l -{ d ^f- ^Q^) 2 , (2.1) 

where A corresponds to the scale at which the interaction term becomes strongly coupled. 
This theory arises in the decoupling limit of the DGP model [12-14], as well as in recent 
models of ghost-free massive gravity [15]. Let us assume that matter couples to the scalar 
via an interaction of the form, 

*h (2 - 2) 

where T is the trace of the energy momentum tensor. It follows that, classically, the 
linearised theory around a non-relativistic source of mass M breaks down at the scale, 

/ \ 1/3 

\M~i) A' ^° see ^ S we no ^ e that the field equations can be schematically written 

(2.3) 

where p corresponds to the momentum operator. For a non-relativistic source of mass, 
M, we get the right schematic behaviour by taking 1 T ~ Mp 3 . In the linearised theory, 
p 2 d> ~ — ij— and we obtain 6 ~ jr-p ~ tt~7.- It is clear that the linearised theory breaks 
down when p 2 cp ~ A 3 , and is only valid for r S> ry. Above the Vainshtein scale we can 
explore corrections to the leading order behaviour by expansion in r /ry. 

Gabadadze et al [9] proposed the following classical dual to the cubic galileon theory 

C! = -I (dcf)) 2 - b^d^cj) - \n<f> + A 3/2 V\b^. (2.4) 

A significant feature of the dual formulation is that the strong coupling scale, A, enters 
with positive powers. Thus the dual Lagrangian is well defined even in the limit A — > 0, 



rv 
as 



1 Strictly speaking, for an approximately point-like source, we would take T = —MS^ 3 '(x) and integrate 
up the differential equations. However, our schematic trick yields exactly the same behaviour, as of course 
it had to on dimensional grounds. 
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in contrast to the original Lagrangian. We also note that the two interaction terms are 
related by a Legendre transformation, 

- ^D0(^) 2 -+ -A 3 / 2 V\P, (2.5) 

with du(j) conjugate to —6^ and □</> conjugate to —A. With this observation, it now seems 
natural to include the canonical kinetic term in the transformation. To this end we define, 

■fo-s^— * , *( 1 + £ D *)- (2 ' 6) 

BC 1 

»™ = m> = - F W) ■ <2 - 7) 

The Legendre transform of the Lagrangian is given by 

A 3 . 



/(6 ( 1), &(2)) = KlM + 6 (2) D ^ " £ = ±A3/ V- 6 (2) 6 (1) " T 6 (2)- ( 2 



The ambiguity in the sign reflects the ambiguity in inverting the relations (2.6) and (2.7). 
This stems from the fact that the cubic galileon theory admits two distinct branches. The 
classical dual is now given by 

C! = b^d^ + 6 (2) D0 - bp)), (2.9) 

and describes the same physics. In the presence of the source (2.2), we obtain the following 
field equations 

d^ 1} -nb {2) = ^-, (2.io) 

pi 

d^±A 3 / 2 b . {2)b{1) " =0, (2.11) 
V- 6 (^(i) 
b 2 a 3 

A 3 / 2 — (1) + — = n (2.12) 

These can be solved order by order in an expansion in A 3 / 2 . For a non-relativistic source 
of mass, M, schematically we have 

pb {1) +p\ 2) = r^-p 3 , (2.13) 



^ + A 3 /y|& (2) |-^ = o, (2.14) 
v |£>(i)| 

p 2 <P + A 3 / 2 -^= + A 3 = 0. (2.15) 

'6(2)| 
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It follows that 




(2.16) 
(2.17) 



M p i 



corresponding to a perturbative expansion in ' 



pj \M v l) ~ \rv) ■ 

3 A general theory with self-interactions and its classical dual 

We will now consider general theories involving self- interactions, extending the ideas initi- 
ated in [9] . In the interests of clarity we present our analysis for a general theory involving 
a single field, (ft, of any type (we suppress tensor indices) with interactions all becoming 
strong in the UV, at the same scale A. Our generalization can be applied to theories with 
multiple scales, and multiple fields, and we will sketch how this should be done in sections 
3.3 and 3.4. 

We start with the Lagrangian density for the field (j> propagating on Minkowski space- 
time 

c(d^<p) = c(</>, d^, . . .), (3.i) 

emphasizing once again that we are suppressing tensor indices on the field - eft does not 
have to correspond to a scalar. We will assume that the field is "canonically normalised" 
in some appropriate way, and that the propagator scales like 1/p 2 in the UV. The various 
interactions are characterised by the number of fields involved in the interaction, N, and 
the number of derivatives, D, and will schematically have the form 



d D (j) 



N 



where D + N > 4. This follows from the fact that we only consider theories that are always 
weakly coupled at low energies, so our interactions should disappear in the limit A — > oo. 
We also require that the interactions become subdominant as —t- so we have N > 2. If 
we further assume that the field couples to a source J via an interaction 

then the field equations are schematically given by 

^ + E^V'- 1 ^. (3.2) 
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Now we can probe the low energy physics by taking A — > oo, in which case we have a 
good linearised theory with p 2 (p ~ J. To probe the high energy physics, we must take the 
opposite limit A — > 0. Which term dominates the dynamics? Naively one might expect the 
dynamics to be dominated by the interaction containing the largest power of 1/A. However, 
things are a little more subtle than that. It turns out that the dynamics is dominated by 
the term (or terms) with largest t, where 

D+N-A , , 

To see this, suppose that the jth interaction dominates the dynamics at high energies. It 
follows that as A —7- 0, <j> oc A' J , since J is independent of A. At the level of the equations 
of motion, the fcth interaction now scales as 



k t 3 \ N k -1 



Since tj > tf~ for all k, this will not diverge as A — > 0. 

Now, if the largest value of t occurs at finite values of N, then it is possible to identify 
the dominant UV behaviour and to expand around it. This controls whether or not we can 
find a classical dual that admits a useful perturbative description. Therefore, assuming 
that 

t = tmax for finite values of N only (3-4) 
we proceed to dualize the theory. 

3.1 The Legendre dual 

To dualize our general theory, we simply compute the Legendre transform for the La- 
grangian. To this end we define, 

a ( o) = W (3-5) 
dC 

a m = dd^> (3 - 6) 

°S = dd^4' (3 ' 7) 



The Legendre transform of the Lagrangian is given by 

/(«(*;)) = O(o)0 + a^d^cj) + a^d^cf) + ...-£. (3.8) 

The precise form of this depends on the inversion of the relation a(&) = q'qI^ - This may 
be multivalued, as it was for the cubic galileon. In any event, choosing some particular 
branch for the inverse, the dual theory is given by 

£' = a (O )0 + a^d^cf) + a^d^d^ + ... - /(a (fc) ), (3.9) 
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with the following field equations 

O(o) - 9/iflJi) + dpdva'Z + ... = - J, (3.10) 

df (3.11) 



da {0) 



W = ^< ( 3 - 12 ) 
ua {\) 

df 

d^d u cf) = — - w , (3.13) 

aa (2) 



Now the original Lagrangian C (d^<f) is such that all interactions remain weakly coupled 
at low energies, in the limit A — > oo. It turns out the dual Lagrangian is well behaved in 
the opposite limit, A — > 0. To see why, we note that schematically 



1 / A \ N i~ 1 



This is inverted to find that the dominant A scaling is (f) ~ A tmax , from which it follows that 
d^(f> ~ A tmax O( k \a) where the operators remain well behaved in the limit A — > 0. 
The Legendre transform is therefore given by 



f(a) ~ A tmax 



A 



£( A W (A:)( a )) 



(3.14) 



Now C(d^(f>) ~(j)J2iP D * (wtY 1 1 and so 



i 

This is well behaved as A — > 0. It follows that the Legendre transform is well behaved as 
A — > 0, provided t max > 0. This is indeed the case since we know that the quadratic term 
has t = 0. We therefore confirm our assertion that the dual theory is well behaved in the 
limit A -> 0. 

We see that the Legendre transform causes the expansion parameter to be inverted, 
the original theory working best at large A, with the dual working best at small A. There 
is a characteristic scale depending on J and A that acts as a pivot about which the duality 
is performed. This is precisely what you mean by the Vainshtein scale in certain modified 
gravity scenarios. On one side of the pivot we have the standard linearised theory with 
corrections that go like negative powers of A. On the other side we have the leading order 
short distance dynamics with corrections going like positive powers of A. The dual theory 
gives us the means to study the latter using ordinary perturbative methods. 

The dual theory describes exactly the same physics as the original theory. This is 
manifestly true when one is considering excitations due to a source. When one is interested 
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in freely propagating modes care must be taken to perform a non-covariant decomposition 
of the conjugate variables 2 , as emphasized in [9]. 



3.2 A practical method for finding the dual 

Now, in general, one cannot explicitly invert the relation = g q^)^ - For example, this 
is already true for the full galileon theory, including quartic and/or quintic interactions. 
[11]. This makes it difficult to find the dual theory using the Legendre transform method. 
Fortunately, however, we may use Lagrange multipliers to arrive at an equivalent dual 
theory with the same useful properties; we will now describe that method. 

Introducing some Lagrange multipliers, C(k)i an d auxiliary fields, An.}, we begin with 
a new Lagrangian, 

C" = C( O )(0 - Ao)) + Cf 1} (<V " + Cf 2 )(d/A</> - A {2)IMU ) + ...+ £(A (fc) ).(3.15) 

This is obviously entirely equivalent to our starting Lagrangian, C, given by (3.1). However, 
the equations of motion for A^) now correspond to constraints that we can use to integrate 
out the Lagrange multipliers. In particular, we find 

c "» = W^' (316) 

c « - 4b (3 - I7) 



£" = ^<t> + 7^-9^ + W + • • • (3-19) 



Plugging this back into the action we obtain, 
dC , dC n , dC 

, . . . dC dC dC 

+C{A {k) ) - A {0) — - Ai),^ " A ^8A^ u - ■ ■ 

This is now of the same form as (3.9), but by using the variables An.} instead of an,} we 
are able to get an explicit expression for the dual Lagrangian. Note that 

dC 

a (t} = gj-, (3.20) 

which is difficult to invert in general. 

As it stands, there is no guarantee that the Lagrangian C" is well behaved as A — > 0. 
We can fix this by rescaling the auxiliary fields An.} = A tmax A^y To see why this helps 
consider the generic interaction 



We thank Gregory Gabadadze for pointing this out. 
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where D = ^ irii and N = Y2% n i ■ It follows that 

where we recall that i = D % N ^ 4: ■ Rescaling our variables, we obtain 

AW (— ) rK) c ^ (— J 3^ n c 

Since t max > i and i max > 0, it is clear that as long as we replace A/ k \ with An^, then both 
£(A( k j) and are wen behaved as A — > 0. As a result, C" is also well behaved in this 

limit. 

Let us see how this method works when applied to the cubic galileon. Note that the 
canonical kinetic term has t = 0, while the interaction has t = 3/2, so we have t max = |. 
Before rescaling we have 

= (-iU^fD - a 1d) v - js A h^ + \ A h + < 3 - 21 ) 

If we set ^4 = h?/ 2 A we obtain the dual Lagrangian 

C" = (-2i (2) if 1) - A 3 / 2 i^) 3 M - 4 )D0 + 1 A 3 i^ + 2A 3 / 2 i (2) i 2 1} , (3.22) 
which is clearly well behaved as A -> 0, as desired. 
3.3 Multiple scales 

Our method for dualising can be easily adapted to deal with theories with more than one 
strong coupling scale. To illustrate how, we consider a theory that depends on two scales 
A <C A, which we may schematically write as, 



3 



We see that we have explicitly separated the interactions into two families - those that 
become strong at A, and those that become strong at A. Again, in order to guarantee 
that all interactions remain weakly coupled at low energies we assume D + N > 4, and 
D + N > 4. We also assume N, N > 2 in order to ensure that the interactions become 
subdominant as <j) — > 0. 

For a source interaction of the form <j)J, the equations of motion are schematically 
given by 

f 2 * + E mk^^- 1 + +E jB^ry-f"'- 1 ~ < 3 - 24 > 

* 3 

To probe the low energy physics we simply take A, A — > oo, and truncate to the linearised 
theory, p 2 ^) ~ J. We can perturb about the leading order solution using inverse powers 
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of A and A. At higher energies we now have two distinct regimes. We first encounter an 
intermediate regime by taking A — > 0, and A — > oo. In contrast, the high energy regime 
is obtained by taking both A, A — > 0. The importance of multiple scales and multiple 
classical regimes was emphasized in general galileon theories in [16] 

Let us first consider the intermediate regime. We can obtain corrections to the leading 
order behaviour in terms of positive powers of A and negative powers of A. Introducing both 
t = g jv_i~ 4 an< ^ * = ^^--i" 4 ' ^ * s c ^ ear fr° m our previous discussion that the intermediate 
scale dynamics is dominated by the term (or terms) with largest t (with the i playing no 
role). To see this explicitly note that the leading order behaviour has 4> ~ A tmax . At the 
level of the field equations, the other interactions go like 

/ \ t \ N-l / ,t \ N-l 

I A <- rnax \ _ / A i-max x 



IT J ' p V a* 

Now since t, i > 0, and t max > t, it is clear that none of these terms diverge as A — > 
0, A oo. 

It is hopefully now obvious how to obtain a suitable dual in the intermediate regime. 
The key point is that the interactions 

1 d D 6* 



f^D+N-A 



are already in the correct form to admit an expansion in terms of inverse powers of A. 
Thus we leave these interactions alone, and focus on dualizing the truncated Lagrangian 

truncated ~ ~ W) 2 + £ — ^ 0°*^. (3.25) 

i 

This can be obtained using either Legendre transforms, or our "practical" method, with 
tmax playing the same critical role as outlined in previous sections. We denote the result- 
ing "truncated' dual as C,' tTuncate( i and note that it is well behaved in the limit A — > 0. 
Combining it with the other interactions, we obtain the following dual for the full theory 

^ = ^truncated + ^- — ^-truncated-, (3.26) 

~ truncated + £ p^^" ^ 
j 

This theory is well behaved as A — > and A — > oo, and beyond this limit we can expand our 
classical solution in positive powers of A and negative powers of A. In a modified gravity 
model with two Vainshtein radii [16], ry and ry < ry, this expansion will ultimately be 
equivalent to an expansion in r/ry, and fy/r. The expansion works well in the intermediate 
regime ry < r < ry. 

We now turn to the high energy regime. Corrections to the leading order behaviour 
are now expressed in terms of positive powers of both A and A. Which term (or terms) 
dominates the dynamics? To answer this we must introduce i max , the largest of the t's, in 
addition to t max , the largest of the t's. The term that dominates the high energy regime 
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tmax qq then 

. In contrast, if 



^34, with 



To see why such terms dominate when they do, we consider the two possibilities sepa- 
rately. If A tmax j A tmax is not divergent, then the A-type interactions with t = t max dominate 
and we have (p ~ A max . To verify that this is correct, we need to show that none of the 
other interactions lead to divergences in the field equations in the desired limit. Indeed, at 
the level of the field equations, the interactions go like 



Since A tmax / A tmax -» oo and i max > t , it is clear that this will not diverge as A, A — > . 

Now consider the alternative scenario, in which, A trnax /A t,nax is not divergent. Using 
an entirely analogous argument one can easily prove that the A-type interactions with 
i = i max dominate. 

It should now be clear how to take arrive at the dual theory in the high energy regime. 
One can simply take the Legendre transform of the full Lagrangian. The dominant terms 
in the expansion of the transform then depend on the ratios of A tmax and A t,nax in the limit. 
In any event, the resulting theory is well behaved as A, A — > 0, and admits a perturbative 
expansion in positive powers of A and A. 

In applying our "practical" method, we need to be sure to rescale the Aha appropri- 
ately, depending on the dominant interaction. In particular, if A tmax /A tmax is not diver- 
gent, we introduce Ae^ = A tmax A^, whilst if A tmax / A tmax is not divergent, we introduce 
Ag^ = A tmax An i y In both cases, the resulting dual theory is guaranteed to be well behaved 
as A, A 0, and to admit a perturbative expansion in positive powers of A and A. In a 
modified gravity model with two Vainshtein radii, ry and fy < ry, the dual theory would 
lend itself to an expansion in both r fry and r/fy. 

The generalization of the ideas presented in this section to theories with even more 
scales should now be obvious. 

3.4 Multiple fields 

We shall now explain how our method should be generalized to deal with more than one 
field. We will assume a single strong coupling scale for brevity, so that the theory may be 
schematically written as 





max — 



> t. The second term 





(3.28) 



a 



- 11 - 



where the index a labels the field and Di = Yl a ^h°" ^ = Ea^V' We have that 
Di + Ni > 4, as well as iVj Q > and Ni > 2. The first condition is required for the theory 
to remains weakly coupled at low energies, while the two latter conditions are required 
in order to guarantee that the interactions become subdominant as <j) a — > 0. We assume 
interactions of the form <f> a J a , so that the equations of motion can be schematically written 
as 



pfy + Z^-zr 1^+^-4 ~ J /3- ( 3 - 29 ) 

i ' 

Now, generically, if a particular interaction (or interactions) dominates the dynamics as 
A 0, then we expect <p a to scale the same way for each value of a 3 . We then claim that 
generically the term with largest ij = — ^^-^ dominates the dynamics. To prove this 
we must show that none of the other interactions will give a divergent contribution to the 
equations of motion as A — > 0. To this end, we note that if our claim is true, each field 
scales as (f) a ~ A tmax . At the level of the equations of motion, the interactions will now go 
like 

\t max ^Di+Ni-A P \ ft* J ' [6.6U) 

where we have used the fact that Di = X)a-^i,<*> = Yla-^ha- Since t max > U and 
Ni > 1 it is clear that this does not diverge as A — > 0. 

It is now clear that generically we can take the dual of this theory in complete analogy 
with the single field case. Again, t max plays a critical role, particularly when applying the 
"practical" method. 



4 The full galileon theory and its classical dual 

As an example of our method, we now consider the full galileon theory 

n=4 

^ = E TS^fr -tf^ 1 ^ • • • ^d^, (4.1) 

n=l (n) 

where each a n = 0(1), and in principle we have a hierarchy of as many as three different 
scales, A( n ), n = 2,3,4. The interaction terms have t n = which will be important in 
establishing how to rescale our conjugate variables. 

Our aim is to find the dual theory that is well behaved as A( n ) — > 0. Adopting the 
practical method, we first arrive at an equivalent Lagrangian, 

3 This is easily seen by taking the ratio of the (3 = /3i and f3 = fii equations of motion. Assuming the 
term with i = j is dominant, we have w 3 '^ 1 x^ 2 - ~ -r^- The right hand side is independent of A, and so 
neglecting the special case where some of the Nj^ vanish, we conclude that <j>p 1 ~ 0,32 
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C' = V -^Mft ...S^A,^ . . . A,^, n (4.2) 

/ j a 3n— 3 y Mi (2)i/i (2)y n v / 

n (n) 

+ E ^3^(0)4 1 ---CJ( 9Ml ^^)^(2 W )t 2 • • • A fa 

n (n) 

/ > ^3n— 3 W Ml (2)fi (2)i/ n 

n (n) 

We need to rescale the variables according to the rules outlined in section 3.3. How 
we do this depends on the ratios of the following in the limit where A( n ) — > 0, 

A*, A*, A*, (4.3) 

where we recall that t n = 3^p. Let us assume, that we are interested in the case where 
A/""\/A/™s is not divergent in this limit for some particular choice of ra* and for any n. 

Then we perform the rescaling A = ^A 

* = (ijg) *fa - *fa (4-4) 

V l \n) J 

/A tn * \ n 

+£ n M # i (O )*...% ] (^^0)i (2 %...i (2 t 



>) 



V A (n) / 



^1 X^nJ J Cl /l Mn 



Because we are assuming that A,^*s/A^\ does not diverge as we take the A( n ) — > 0, it 
follows that this dual action is well behaved in the limit. 

Assuming that matter couples as in the cubic galileon case, we have the following 
equations of motion in the dual theory, 

n=4 

T,^^..^ a fa... a fa 



n=l 
n=4 



+ Y,na n {\ [n) Y5^..^d^d Ul [A [0) Afa...Afa) + J = 0, (4.5) 

n=l 
n=4 

E«"( A w) n *-CJ {^^-KnA^) ^••• i (2K=°' ^ 
n=l 

n=4 

E ^(\n))^fe...tfi (2 ^ 2 • • • A fa 
n=l 
n=4 

+ J>(»- l)a B (A (B) )M ( o ) «ft...%l^M^, ••• A fa 



n=l 
n=4 



£ n 2 a n A;; :) (A n )M (0) 5fe...^Ji (2 ^ . . . = 0, (4.7) 

n=l 
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where A( n ) = A/^\/A^\ . We now find the background solution in the following high energy 



(«*)' («) 
limit: A( n ) — > with A( n ) 
the background becomes 



1 n = n* 
n / n* 



. Taking the source to be J ~ w~P 3 ^ as before, 



1(2) 



l/n* 



-4(0) ~ A (2)P~ 



M P i, 
M N 1/n * 



P 



9 f /MX"''"* 3-2n* 



K n \ A (2)P~ ~ A "\ — p~ 

K) W (n*) ^ Mp i / F 



1/n* 



(4.8) 
(4.9) 
(4.10) 



with the higher order terms being by expansions in A/™*% and in A'™^ for re 7^ n*. These 



correspond to expansions in 



and 



3(n-n*) / M 



M 



M 



-in* /3 



3(n* — n) 



respectively. The theory admits up to three critical radii given by 



n / re* 



M \ \ 3(n-n») 



\n) 



re = re* 



re 7^ re* 



thereby generalising the Vainshtein radius for the multiscale theory, as expected [16]. 



5 Discussion 



Whenever classical perturbation theory breaks down at some particular scale, to continue 
making predictions beyond that scale one would like to have a dual theory. This should 
describe the same physics, but admit a perturbative description that works best in the 
opposite regime to that in the original theory. This is exactly what Gabadadze et al [9] 
achieved by identifying the duals to two particular classical theories that exhibit Vainhstein 
screening. In each case, perturbation theory works best above the Vainshtein radius in the 
original theory, and below the Vainshtein radius in the dual. 

In this paper, we have recognised these examples as being nothing more than the 
Legendre transform of the original Lagrangian. This has enabled us to generalize the idea, 
and outline how one can find the classical dual for a much larger class of self-interacting 
theories. We have also presented a more "user- friendly" method for finding the dual for 
the case where it is difficult to compute the Legendre transform explicitly. 

Of course, a dual is only of any use if it admits a complementary perturbative descrip- 
tion, as in the examples given in [9]. If the original theory admits a good perturbative 
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description as some dimensionful scale A — > oo, this amounts to the dual theory admitting 
a good perturbative description in the opposite limit A — > 0. That this should happen 
places certain restrictions on the form of the theory one can successfully dualize. The key 
ingredient is that there must exist a finite interaction (or interactions) that dominates the 
classical dynamics as A — > 0, and about which we turn the theory on its head. This logic 
applies even when there are multiple scales or multiple fields. 

We have presented the classical dual of the full galileon theory as an example. Of 
course, there are many more examples one could consider, and one can find their duals 
using the methods we have discussed. We have also shown explicitly that a useful dual can 
only be found if certain specific criteria are met: namely that there exists a finite Appoint 
interaction (or interactions) with largest t, where t = D ~jy^ 4: and D is the number of 
derivatives. 

Let us end with an example of a Lagrangian for which we cannot find a useful classical 
dual. Consider the Einstein-Hilbert action, written in terms of the metric expanded around 
Minkowski space, g^ v = rj^ u + j^—h^. Schematically we have 

pi 

/oo ^2 n 
i=3 M PI 

As is well known, this theory is well behaved as M p i — > oo, as the interactions vanish. Each 
interaction has Di = 2 derivatives, Ni = i fields and so t{ = g '^.^ 1 ~ 4 = — f^f ■ Recall that 
the dynamics as M p i — > is dominated by the term (or terms) with largest tj. However, 
this does not occur at a finite value of i, and so we cannot find a classical dual that is well 
behaved as M p i — > 0. This does not mean that one cannot find a classical dual to GR. 
Such a dual can be found but only if we describe the dynamics using something other than 
h^y. This is currently a work in progress [17]. 
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